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1.  INTRODUCriON 


The  DYNASPHERE  computer  code  for  the  solution  of  system-generated  elec- 
tromagnetic pulse  (SGEMP)  problems  in  spherical  geometry  is  described  in 
this  report.  Emphasis  is  given  to  the  practical  details  necessary  for  the 
effective  operation  of  the  code.  A brief  description  of  the  code  is  out- 
lined, and  then  specifics  such  as  computer  requirements  essential  to  its 
operation  are  discussed.  Input  card  requirements  are  detailed,  a sample 
problem  is  treated,  and  the  code  output  interpreted.  Many  prominent  code 
variables  are  defined  in  a glossary.  All  terms  necessary  to  the  effective 
application  of  DYNASPHERE  to  typical  SGEMP  calculations  arc  defined.  De- 
tailed descriptions  of  physics  and  modeling,  as  well  as  ci>de  checkouts,  are 
found  in  Reference  1. 


1,  T.  N.  Delmer  et  al.,  "SGEMP  Phenomenology  and  Computer  Code  Devel- 
opment," DNA  3653F,  November  11,  1974. 
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PRSCEDINO  PAOS  BLMfC.NOT  flLMSD 


2.  t)YNASPHi:Ri:  ohscription 


DYNASMPFRl-  treats  the  SCILMP  two-dimensional  problem  for  axisymmelric* 
electron  emission  from  concentric  perfect  ly-condiict  inn  spheres.  Ilmission 
electron  enerRV  spectrum,  pulse  shape,  and  spatial  distributions  must  be 
specified  to  the  code.  Maxwell's  equations  are  solved  to  obtain  fields 
which  act  on  the  electrons.  The  electron  motion  is  treated  by  following 
particles  of  charge  through  the  s|>atial  grid.  The  calculation  is  self- 
consistent  in  that  fields  modify  electron  trajectories  which,  in  turn, 
modify  the  fields. 

Electron  emission  can  occur  from  the  inner  or  outer  spheres,  although 
presently  the  code  will  not  emit  from  both  spheres  simultaneously.  Energy 
spectra  and  spatial  and  angular  distributions  are  arbitrary.  The  emission 
current  pulse  is  specified  with  a trapezoidal  shape,  with  arbitrary  rise 
and  fall  times  and  pulse  length.  No  t ime-retardat ion  of  the  emission  cur- 
rent is  modeled.  This  can  he  a considerable  deficiency  for  small  ratios 
of  pulse  length  to  light  t imes-of- fl i ght  across  object  dimensions,  or  in 
calculations  where  surface  currents  within  the  emission  region  are  required 
(Ref.  2 j . Worst-case  surface  currents,  generally  located  just  outside  the 
emission  region,  occur  in  situations  where  moderately  long-duration  pulses 
are  coupled  with  simultaneous  emission.  Minor  code  modifications  could  be 
made  to  retard  electron  emission  for  very  short-pulse  problems. 

Particles  of  charge  representing  large  numbers  of  electrons  arc  used 
to  represent  photo-electric  emission  from  surfaces.  These  particles  can 
be  injected  with  discrete  or  random  energy  and  angular  distributions. 
Electron  emission  at  discrete  spatial  positions  are  required  presently. 

The  normal  configuration  used  in  the  calculations  consits  of  concen- 
tric cylinders,  although  an  isolated  sphere  can  be  treated  in  the 

2.  E.  P.  Wennas,  S.  Rogers,  and  A.  .1.  Woods,  "Sensitivity  of  St'.l'MP 
Response  to  Input  Parameters,"  IEEE  Trans.  Nucl . Sci.  NS-22 , No.  0, 

December  1975. 

•* 

The  axis-of-emission  symmetry  is  defined  by  the  direction  of  photon 
propagat i on . 
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qiias  i-stat  ic  section  of  the  code.  OYNASPUrRII  performs  the  calculation 
for  the  electric  and  magnetic  fields  produced  by  the  electrons  using 
either  the  full  Maxwell's  equations  or  the  quasi-static  Green's  function 
approach,  depending  on  input  options  specified.  The  latter  a]>proach  was 
employed  in  the  forerunner  code  r.sl’Hi.Ki;  (Kef  1)  due  to  its  ease  of  imple- 
ment.at  ion.  The  method  is  valid  under  conditions  of  long  pulse  rise  times 
and  low-energy  electron  spectra  (Ref  .').  It  has  been  retained  as  a subset 
of  the  DYNASPIlliRi;  code,  but  is  not  discussed  in  detail  here  because  of  the 
former  documentation.  The  full  Maxwell's  equation  treatment  permits  cal- 
culation of  la*:e-timc  current  ringing  in  the  structure  and  provides  general- 
ly more  accurate  results. 

Spherical  coordinates  are  employed  for  the  spatial  zoning  in  liYN/VSI’IILRh. 
Excellent  resolution  of  the  region  surrounding  spherical  objects  can  be 
obtained  conveniently  with  this  coordinate  system.  Spatial  zoning  can  be 
variable  in  the  radial  direction  hut  must  be  azimuthal ly  symmetric.  Func- 
tions to  specify  radial  zoning  to  the  code  are  employed  internally.  Con- 
stant or  variable  zoning  may  be  requested.  In  the  latter  case,  sizes  can 
vary  only  mildly  (typically  a factor  of  2 increase  from  the  inner  sphere  to 
the  outer  sphere)  or  rapidly  (typically  a factor  of  15  increase  from  inner 
to  outer  sjihere) . Strong  variations  are  generally  employed  for  high  space- 
charge-  1 imited  (SCI.)  conditions  where  resolution  of  large  field  and  current 
gradients  is  required. 

Forces  acting  on  the  p.irticles  representing  photo-electrons  are  limited 
to  electric  fields  only. 

Two  time  steps  arc  employed  in  dynamic  ca.'.cs.  The  Ma.xwell's  equation 
time  step,  or  "light  time  step,"  is  determined  automatically  by  the  code. 

It  satisfies  the  Courant  stability  criterion  and  is  also  an  integer  divisor 
of  the  particle  or  "electron  time  step.”  Thus,  the  field  equations  may  be 
called  several  times  for  each  particle  position  update.  A minor  code  modi- 
fication is  necessary  to  specify  any  particular  light-time  step. 

Ti  in  PT  Wenaas  and  A.  .J.  Woods,  "Comjiarisons  of  Quasi-Static  and 
Fully  Dynamic  Solutions  for  Electromagnetic  Field  Calculations  in  a 
Cylindrical  Cavity,"  IEEE  Trans.  Nucl.  Sci.  NS-21 , No.  b,  December  1974. 
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A convenient  summary  of  the  DVNASPIIliRH  culculat  ional  sequence  is 
found  in  the  flow  chart  in  Figure  1.  Notice  that  elect ron- emission 
information  is  obtained  from  a separate  source.  Also,  notice  the  optional 
electron-trajectory  "movie*’  capability.  These  movies  can  be  obtained  by 
specifying  an  input  option  and  saving  the  resulting  particle  information 
on  tape  for  later  treatment  by  the  MOVIF.  code  (Ref.  1). 


4.  A.  J.  Woods,  T.  N.  Delmcr,  and  M.  A.  Chipman,  "The  Arbitrary  Body 
of  Revolution  Code  (ABORC)  for  SC.HMP/lF.Mr,"  INTF.L-RT  8141-028,  April  197b. 


6 


RT-13703A 


Figure  1.  OYNASPUHRF.  flow  chart 


3,  llYNA-irHURK  COMPUTI  R HI  qUIRl>«LNl  S 


DYNASPHJiRH  is  a l•OR^R,■\^- IV  computer  proRrara  of  al>out  4000  cards, 
including  both  the  dvTiamic  and  ipiasi -stat  ic  segments  of  the  code.  No 
machine  language  coding  is  employeil.  The  code  operates  on  the  ClHI  7000 
computer.  Core  reijuirements  are  12000^^^  words  of  small  core  .and  IOO.OOOjq 
words  of  large  core.  Three  fast -access  tiles  ;irc  also  rei|utrcd  during 

execution,  and  four  arc  optional,  depending  on  code  input  options  chosen.  J 

The  RUN  compiler  is  currently  used,  although  conversion  to  the  faster  i 

FTN  compiler  should  be  straightforward  at  this  time.  Run  times  vary  from 
1 to  30  minutes  of  central  processor  time,  depending  on  problem  condi- 
tions. Core  must  be  preset  to  zero. 
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4,  DliSCKlPTION  01  nil  IM'inS 


Petailcii  descrijit  ions  of  input  >juaiu  i»  Jos  lequired  by  KYNASI’llbRi;  are 
listed  in  this  section.  \ariablt“  names  and  their  physical  or  calcula- 
tional  si^ini f icance  are  given,  along  with  f»r  .its  for  reading  them  into 
the  code. 

A brief  list  of  definitions  peculiar  to  1>V\ASI’III.RI;  is  given  in  Table 
1.  This  information  expedites  communication  «f  detailed  cotie  quantities. 
Variable  limits  are  found  in  Table  2;  these  help  the  user  to  stay  within 
the  numerical  and  physical  limitations  of  the  code.  The  positions  on 
the  grid  where  quantities  are  evaluated  are  given  in  Table 

In  addition  to  the  detailed  input  descriptions  listed  here,  the  user 
will  find  a complete,  abbreviated  input  description  in  the  code  listing 
itself.  The  descriptions  are  at  the  beginning  of  the  main  prognun  (called 
MAIN).  Once  the  progr.immer  understands  the  basic  functions  of  the  inputs, 
the  abbreviated  manual  is  probably  the  most  convenient  to  use.  Also,  cer- 
tain variables  are  not  defined  in  the  input  descriptions  in  this  report 
because  they  pertain  either  to  debugging  the  older  quasi-static  version 
of  the  code  or  to  optional  editing  feature’s.  The  abbreviated  manual  con- 
tains definitions  of  these  quantities.  These  items  can  be  ignored  in 
standard  SGF.MP  calculations. 

Input  card  descriptions  appear  at  the  end  of  this  section.  The  card 
numbers  shown  actually  mean  "card  type".  If  an  array  requires  more  numbers 
than  fit  on  one  card,  it  is  continued  on  the  next  card.  .An  array,  as 
opposed  to  a single  variable,  is  identified  quickly  from  the  format  infor- 
mation. A number  appearing  in  front  of  the  format  (A,  R,  or  I)  indi- 

cates an  array,  while  lack  of  a number  indicates  a single  .alue. 

Variables  may  have  more  than  one  name.  All  names  are  given  for  con- 
venience. The  words  "FDIT"  and  "nEBUG"  appearing  at  the  beginning  of 
descriptions  imply  that  the  variable  is  used  either  for  editing  or  for 
debugging  purposes.  These  titles  help  the  user  to  scan  the  descriptions 


more  rapidly.  The  designations  "MAX  and  "MAX  SO"  give  the  maximum  value 
and  the  maximum  number  of  values  a variable  can  have,  respectively.  Min 
imum  numbers  of  values  are  given  by  "MIN  .SO"  wliere  appropriate.  Default 
values,  set  by  the  code  when  no  value  is  read  in,  are  also  noted  where 
appl i cable. 


Table  1 

DEFINITIONS  OF  TFR.MS  COMMONLY  USED  IN  THIS  MA.SUAL 
.-VND  IN  DY.SASPlIERi:  OUTPUT 


Particles 


Mini-print 


2-D  print 


Dynamic  calculation 


Quasi-static  calculation 


Large  numbers  of  photo-electrons 
grouped  into  a single  charged 
pa rt icle. 

Short  printout  (3  lines)  giving 
summary  quantities  helpful  in  ana- 
lyzing the  calculation.  Can  be 
printed  out  at  times  independent 
of  the  large  2-D  prints. 

Large  printout  giving  spatial 
distributions  of  fields  and  cur- 
rents, as  well  as  other  pertinent 
information.  Can  be  printed  out 
less  often  than  the  small  mini- 
prints. 

Indicates  full  Maxwell's  equations 
solution  for  fields,  as  opposed  to 
the  quasi -static  approximation. 

Indicates  fields  calculated  from 
Maxwell's  equations  with  the  time- 
derivative  terms  removed.  Green's 
function  technique  is  used. 
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MINIMUM  AND  MAXIMUM  VALUES  OF  DYNASPliERE  VARIABLES 


Min.  N'o. 

Max . No . 

Input  Variables 

Zone  number 

1 

3 j 

50  radial  I , 

1 / dynamic 

20  angular^ 

20  radial  ) quasi- 
20  angular/  static 

Time  steps 

- 

- 

Emission  zones 

1 

50 

Energy  distributions 

1 

50^* 

Energy  bins 

1 

50 

Time  histories 

1 

1 

Angular  distributions 

1 

50 

Angular  distribution  bins 

50 

2-D  prints 

- 

- 

Mini-prints 

- 

- 

Calculational  Variables 

Total  number  of  particles 
emitted  in  a given  time  step 

0 

2000 

Number  of  particles  being 
followed  in  a given  time  step 

0 

20000 

^Only  1 angular  and  energy  distribution  per  emission  zone  are  permitted, 
however. 


Tabic  3 


Ki;i,ATl\T  POSITIONS  ON  OKI!)  MIliRP  l•■ll;l.l'S,  QIRRIiNTS, 
AND  CHAlU'.l'.S  ARP.  CAl.CDl.ATi:!) 

Positions  c and  b indicate 
zone  center  and  boundary, 
respectively.  Tlicsc  posi- 
tions pertain  to  dynamic 
calculations  only. 


A.\  i a 1 

Radial 

Quant i ty 

Posit  ion 

Posit  ion 

r. 

r 

c 

b 

P 

u 

c 

c 

c 

J 

r 

c 

b 

1) 

c 

Charge  in  zone 

b 

b 

Surface  current 

c 

0 

(a)  0 is  the  polar  angl 

c relative  to  the  photon  pro]iagation  direction 

(b)  (J)  is  the  azimuthal 

angle  relative  to  the  plioton 

propagat ion  d i rect ion 

INPUT  CARDS 


Card 

Number 

1 

2 

2 

2 

2 

2 

2 


Columns  Variable 

(format)  Name(s)  Description 


1-78 

(13A6) 

T1TU-. 

Comment  card 

1-2 

(12) 

lOPT(l) 

(MEITUOD) 

Quasi-static  calculations  only:  gives 
Green 's  function  soui’cc  chai'ge  type 

1 Source  cliargcs  are  rings 

2 Source  charges  smeared  out  over 
sjiat  ial  zone 

3-4 

(12) 

IOPT(2) 

(IRPINT) 

PDIT:  Limits  printout 

-1  Minimum  jirintout  (no  particle 
emission  cliaractcr i st ics) 

0 Particle  emission  character- 
istics 

1 Green's  function  print  if 
cpiasi-stat  ic 

S-b 

(12) 

IOPT(3) 

(INPUTR) 

Determines  emission  particle  initial 
radial  positions  in  meters 

0 limit  ])articles  from  radial 
pos i t ion 

r = 0.999* (RSAT  - ri/2)l,  dynamic 

r = 0.999*RSAr,  quasi-static 
>0  limit  jiarticles  from  radial  posi- 
tion INPUTR  (MAX  S outer  sphere 
radius) 

7-8 

(12) 

IOPT(4) 

(ITDIM) 

Do  not  e.xecute  problem  of  >0. 

Allows  for  stopping  after  Green's 
function  calculation  in  quasi-static 
cases.  Code  does  not  read  any  infor- 
mation beyond  the  zoning  cards  (cards 
.3-8).  Pennits  checkout  of  Green's 
function  on  various  analytic  charge 
densities.  See  10PT(20). 

9-10 

(12) 

10PT(5) 

(TDUMP) 

DEBUG;  Error  off,  to  get  core  dump 
at  end  of  jirohlem  if  >0. 

11-12 

(12) 

lOPT(b) 

(IGRID) 

Quasi-static  only; 

0 Source  and  field  point  grids  are 
the  same 

1 Source  and  field  point  grids  will 
be  different 

Card 

Number 

Columns 

(format) 

Variable 

Name(s) 

Descript  ion 

-y 

15-14 

I0PT(7) 

DEBUG:  Permits  Green's  function 

(12) 

(IQS AT) 

checkout  with  no  contribution  from 
net  charge  on  inner  sphere 

2 

15-16 

(12) 

I0PT(8) 

DEBUG:  Quasi -static  only. 

A number  of  test  problems  to  run 

checking  out  Green's  function. 

If  >1,  code  reads  in  I0PT(8)-1  new 
values  for  the  variables  IOPT(7,12, 
13,14,20)  in  sequence  after  zoning 
cards  (.cards  5-8).  Permits  redefi- 
nition of  various  analytical  charge 
densities.  See  IOPT(.20)  and  also 
abbreviated  user's  manual  in  MAIN 
program . 


2 

43-44 

(12) 

10PT(22) 

(NESKIP) 

Emit  particles  every  NESKIP+1  time 
steps.  Default  = 0. 

2 

47-48 

(12) 

I OPT (24) 
(NTSKIP) 

Print  out  2-U  prints  every  IOPT(24) 
particle  time  stci)s.  Overridden  by 
DTPRNT  (card  16). 

2 

40-50 

(12) 

IOPT(24) 

(I  EMIT) 

0 Define  emission  i)article  angu- 
lar bins  relative  to  axis  of 
coordinate  system. 

1 Define  emission  particle  angu- 
lar biiKs  relative  to  surface 
normal . 

2 

51-52 

(12) 

10PT(26) 

(I0P26) 

(NEPTS) 

Number  of  electron  emission  zones. 
If  zero,  emission  zones  same  as 
angular  zones  defined  by  card  6, 
and  cards  11  and  12  need  not  be 
read  i n . MAX  = 50 . 

2 

53-54 

(12) 

IOPT(27) 

(NTH) 

Number  of  emission  electron  angular 
lar  bins.  MAX  = 50. 

2 

55-56 

(12) 

I0PT(28) 

(NSPD) 

Number  of  emission  electron  energy 
bins.  MAX  = 50. 

2 

57-58 

(12) 

IOPT(29) 

0 Quasi-static  approximation. 

1 Full  .^laxwell's  equation  treat- 

ment . 

-1  Quasi-static  a]iproximation  with 
constant  zone  size  in  both  radial 
and  angular  directions.  Do  not 
read  in  cards  5-8. 


L 
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Card 

Columns 

Variable 

Number 

(format) 

Name(s) 

(iesc  ript  ion 

2 

59-60 

IOPT(30) 

RDIT:  Printout  limiter  in  dynamic 

(12) 

(IPRNT) 

cases  only. 

0 Shnrt  print  - docs  not  print 
out  surface  current  and  charge 
density  on  outer  sphere,  nor 
photo-electron  charge  density 
in  2-1)  prints 

1 Long  print  - prints  out  surface 
currents  and  charge  densities  on 
outer  sphere  as  well  as  photo- 
electron charge  density  in  2-D 
jirint.  Also  prints  out  several 
other  arrays  such  as  charge  in 
each  2one  as  well  as  current 
densities  calculated  in  a man- 
ner different  from  the  Maxwell's 
equation  section. 


2 

61-62 

(12) 

I0PT(31) 

Print  mini-prints  every  10PT(31) 
particle  time  steps.  DPd'AULT  = 1. 

2 

63-64 

(12) 

IOPT(32) 

MOVIF:  option:  Write  particle  posi- 
tion information  to  file  TAPE20 

every  lOP'I'f.'iZ)  particle  time  steps 
if  >0.  If  <0,  plot  particle  posi- 
tions on  jirinter  as  well  as  writing 
file. 


3 

1-2 

KOPT(l) 

Number  of  radial  zones 

(12) 

(NR) 

MIN  = 3 

MAX  = 50 

3 

3-4 

K0PT(2) 

Number  of  angular  zones 

(12) 

(NT) 

MIN  = 3 

MAX  =20 

3 

5-6 

K0PT(3) 

Use  in  quasi-static  cases  only 

(13) 

(NC) 

0 Radial  source  point  zoning  in 
Green's  function  same  as  field 
point  zoning 

>0  Number  of  radial  zones  for 
sources  for  fields  in  Green's 
function.  Sec  IGRID  [10PT(6)] 
MAX  = 20 

3 

7-8 

K0PT(4) 

Like  K0PT(3)  but  for  angular  zones 

(12) 

(NA) 

MAX  = 20 

1.1 


Card 

Number 


Columns 
( format) 


Variable 

Name(s) 


Description 


HDIT:  2-D  print  size  limiter. 

0  Print  out  approximately  10  radial 
zones  in  2-1)  prints  in  addition 
to  special  zones  specified  by 
K0PT(7  and  8) 

>0  Print  out  cv'ery  K0PT(5)  radial 
zones  where  indices  arc  between 
the  values  specified  by  K0PT(7 
and  8).  Di;i-AULT  = 1. 


5 

11-12 

(12) 

K0PT(6) 

EDIT;  Print  out  every  K0PT(6)  angu- 
lar zones. 

DEFAULT  = printout  5 It)  angles. 

3 

13-14 

(12) 

K0PT(7) 

Print  every  radial  zone  up  to  and 
including  zone  number  K0PT(7) 
regardless  of  K0PT(5). 

DEFAULT  = 1. 

3 

15-16 

(12) 

K0PT(8) 

EDIT;  Print  out  every  radial  zone 
with  index  ^ K0PT(8)  regardless  of 
K0PT(5) . 

DEFAULT  = number  of  radial  zones 

Radial  zoninj;  flag  - dynamic  calcu- 
lations only, 

0 Constant  zone  size. 

1 nilBUG:  .Sets  constant  zone  size 
but  uses  I-'DNQR  subroutine.  Per- 
mits checkout  of  that  routine 
when  compared  with  same  calcula- 
tion with  KOI’ 1(01  = 0. 

2 Zone  size  increases  slowly  to 
about  twice  as  much  at  outer 
sphere  as  at  inner  siihcre. 

3 Zone  size  increases  rapidly  to 
about  10  times  as  large  at  outer 
sphere  as  at  inner  siihcre. 

3 19-20  KOPT(IO)  Random  emission  flag; 

(12)  0 |)iscrete  emission. 

>0  Randomize  energy  and  angle  of 
emission  particles  iictwccn  bin 
centers  defined  by  cards  9,  10, 
1,3,  Id,  Reijuires  at  least  2 
bins  for  energy  spectra,  2 for 
angular  distributions.  Different 
values  of  KOPT(IO)  give  different 
random  number  sequences. 


3 17-18  K0PT(9) 

(12) 


3 9-10  K0PT(5) 

(12) 


1( 


Card 

Number 


Columns 

(format) 


Variable 

Name(s) 


Hcser ipt ion 


3 

21-22 

(12) 

KOPT(ll) 

3 

23-24 

(12) 

K0PT(12) 

4 

1-12 

(E12) 

RSAT 

4 

13-24 

(E12) 

B 

5 

1-72 

(6E12) 

Cl 

6 

1-72 

(6E12) 

ALPHA 

7 

1-72 

(6E12) 

RI 

8 

1-72 

(6E12) 

THETA 

9 

1-72 

(6E12) 

ENRG 

10 

1-72 

THETAO 

If  >0.  v»rite  current  densities  in 
radial  and  angular  directions  to 
file  TAPfJl  every  particle  time 
step.  I’crmits  use  by  lumped-element 
modeling  code,  which  generates  time 
histories  of  the  currents. 

Check  for  more  than  1 input  deck  if 

>0. 

Radius  of  inner  sphere  (m) 

Radius  of  outer  sphere.  For  quasi - 
static  calculations  with  single 
sphere  only,  set  RSAT  >10*^  (m) 

Radial  cone  boundaries.  Read  in  for 
quasi -static  cases  only  (m). 

MAX  NO  = 20 

Angular  cone  boundaries.  Read  in  for 
quasi -static  cases  only  (radians). 

MAX  NO  = 20. 

Radial  zone  boundaries  for  source 
charges  in  Green's  function.  Read 
in  in  quasi -static  cases  only  (m) . 

mx  .\0  = 20. 

Angular  zone  boundaries  for  source 
charges  in  Green's  function.  Read 
in  in  (|uasi -stat ic  cases  only 
(radians  ) . 
mx  NO  = 20. 

Initial  energies  of  emission  elec- 
trons. Read  in  10PT(28)  values 
(keV) . 

MAX  NO  = 50. 

TI1F.TA-"ZI:R0" : Initial  directions  of 
emission  electrons  at  each  emission 
point.  Measured  from  surface  nor- 
mal or  z a.xis,  depending  on  IOPT(25). 
MAX  NO  = 50. 


Card 

Number 


Columns 

(format) 


Variable 
Name (s) 


Description 


i 

1 


11 


12 


13 


14 


15 


16 


1-72 

(6E12) 


1-72 

(bF.12) 


1-72 

(6F.12) 


1-72 

(6E12) 


1-72 

(6E12) 


1-12 

(E12) 


nip 


DTH 


ENORMS 


ELN 


FTH 


TMAX 


Angular  positions  of  emission  points 
mcasureil  from  axis  of  the  coordinate 
system.  Read  in  10PT(26)  values  only 
if  IOF’T(26)  X)  (radians). 

MAX  .NO  = .50. 

Width  of  the  emitting  sectors  in  units 
of  cosinetmin)  - cosinc(max)  where 
min  and  max  can  be  thought  of  as 
the  minimum  and  maximum  angular 
positions  of  the  emission  :one  edges. 
Not  a convenient  input.  See  10PT(26) 
for  eas ier-to-use  method.  Read  in 
10I'T(20)  values  if  I OPT  (26)  >0. 
(dimensionless) 

MAX  NO  = 50. 

Normalizing  factor  for  emission  cur- 
rent densit)’  from  each  emission  zone. 
Multiplies  .IPEAK.  Negative  value 
means  use  the  energy  and  angular  dis- 
tributions from  the  previous  emission 
zone  (1  zone  closer  to  the  axis). 

Read  in  10PT(2o)  values  (dimensionless) 

Relative  number  intensity  of  elec- 
trons in  each  energy  group.  Read  in 
as  many  cards  1 1 and  15  in  pairs  as 
there  arc  positive  values  of  FNORMS 
(card  13).  Code  normalizes  spectra 
to  Jl’tiAK  jieak  emi.ssion  current  den- 
sity. Beware  of  these  units. 

( e 1 ec t ron s/ ene rg)-  bin) 

MAX  NO  = 50. 

Relative  number  intensity  of  emission 
electrons  in  each  emission  angular 
distribution  bin.  Sec  also  card  14. 

All  distributions  normalized  by  code 
to  peak  emission  current  density  of 
.JPEAK.  Beware  of  these  units, 
(electrons/angular  bin) 

MAX  NO  = 50 

Maximum  time  of  calculation  (nsec) 
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Card 

Number 

Columns 
( format ) 

Variable 

Name(s) 

Descript  ion 

16 

13-24 

(F.12) 

DF.LT 

Particle  time  step  (nsec) 

16 

25-36 

(F12) 

T1 

Rise  time  of  the  trapezoidal  emission 
current  pulse  (nsec) 

16 

37-48 

(E12) 

T2 

Emission  current  pulse  is  constant 
from  T1  to  T2  (nsec) 

16 

49-60 

(E12) 

T3 

Emission  current  pulse  falls  off 
linearly  to  zero  from  T2  to  T3  (nsec) 

16 

61-72 

(E12) 

DTPRNT 

Print  2-D  prints  every  DTPRNT  nsec 

17 

1-12 

(E12) 

NTOT 

Total  number  of  electrons  emitted  in 
the  emission  current  pulse.  If  zero, 
calculated  from  JPEAK  (card  17). 
(electrons) 

17 

13-24 

(E12) 

JPEAK 

Peak  emission  current  density.  If 
zero,  calculated  from  NTOT  (card  17). 
Either  NTOT  or  JPEAK  must  be  non- 
zero. (amp/m^) 

17 

25-36 

(E12) 

TSTART 

Quasi-static  calculations  only; 

Restart  from  dump  tape  at  time 

5'  TSTART.  See  10PT(17).  (nsec) 

17 

37-48 

(E12) 

DTTAPE 

Quasi-static  calculations  only: 

Write  dump  tape  every  D'lTAPE  nsec. 

See  I0PT(18). 
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5.  nYNASnpiiRi;  S;Vmi'lh  problbi 


A sample  calculation  is  described  in  this  section.  The  physical  prob- 
lem is  described,  and  the  required  input  cards  are  jiiven.  Selected  code 
outputs  are  listed.  Descriptions  of  these  outputs  are  found  in  the  variable 
glossary  of  Section  (>. 

This  sample  problem  is  sufficient  to  illustrate  many  code  features  and 
to  provide  a test  of  DYNASPUHRH  when  it  is  converted  to  a different  computer 
system. 

The  sample  problem  consists  of  a si)here  of  2..i  m radius  enclosed  in  an 
outer  sphere  of  15.2  m radius.  Mono-energetic  electrons  of  energy  7.1  keV 
are  emitted  uniformly  from  one-half  the  inner  sphere  surface.  The  pulse 
shape  is  trapezoidal,  with  45  nsec  rise  and  fall  times  and  55  nsec  full- 
width-at-half-ma.ximum.  The  peak  emission  current  density  is  0.108  amp/m*', 
and  the  electrons  are  assumed  to  be  emitted  radially  outward. 

The  problem  conditions  result  in  a non-space-charge- 1 imited  solution. 
Gradients  will  not  be  steep,  and  the  slowly  varying  radial  zone  size  is 
sufficient  to  give  reasonably  accurate  resolution  of  gradients  near  the 
inner  sphere.  K0PT(9)  is  set  to  2,  and  30  radial  zones  are  used,  result- 
ing in  a minimum  radial  zone  size  of  0.21  m (chosen  by  the  code).  A time 
step  of  2.5  nsec  causes  particles  to  traverse  about  one-half  this  distance 
each  step.  Ten  angular  zones  and  five  emission  zones  are  employed,  with 
the  latter  beiiig  set  by  the  code.  I'he  prolilem  is  run  to  150  nsec,  witli 
large  sinitial  printouts  every  20  nsec.  Particles  are  emitted  every  other 
time  step. 

The  particular  model  chosen  for  this  sample  problem  should  not  be 
regarded  as  a typical  SGHMP  problem.  Simplifying  assumjit ions , such  as 
monoenerget ic  electron  emission  normal  to  the  surface,  may  give  unrealistic 
results,  and  the  practice  of  omitting  particles  every  other  time  step  can 
be  particularly  dangerous  in  siiace-chargo-1 imited  situations.  These  simp- 
lifying assumptions  were  adequate  for  the  purposes  of  this  sample  problem, 
however. 

20 


Input  card  images  describing  the  sample  problem  to  the  code  are 
shown  in  Figure  2.  The  particle  emission  characteristics  resulting  from 
the  electron  emission  specifications  are  shown  in  Figure  3.  iTiese  var- 
iables give  initial  electron  positions,  energies,  and  angles.  A sample 
mini-print  and  a sample  si>atial  print  apjjear  in  Figures  t and  5,  for  a 
time  of  20  nsec.  Limited  information  on  fields,  currents,  particle 
positions,  etc.,  is  given  in  the  mini-print,  while  detailed  spatial 
distributions  are  contained  in  the  larger  print.  .Note  that  only  12  of 
the  30  radial  zones  arc  printed  out  in  the  spatial  distribution  of  fields 
and  currents.  All  30  zones  could  liave  been  obtained  by  setting  K0PT(5) 

= 1. 

A slightly  low  value  for  the  quantity  QIN.IFK  (total  emitted  charge) 
will  be  seen  in  the  mini-print.  The  number  shown  is  about  2S  percent 
lower  than  the  time  integral  of  the  emission  current.  This  is  a conse- 
quence of  choosing  a large  time  step  and  emitting  jiarticles  every  other 
step  (a  procedure  that  is  not  recommended).  By  the  end  of  the  pulse,  the 
emitted  charge  is  much  more  consistent  with  the  expected  value. 
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Figure  2.  OYNASPIiFRI-  sample  problem  input  deck 


Figure  3.  DYNASPHERF.  particle  emission  characteristics  printout  sample 
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Figure  5.  DYNASPIIFRE  spatial  print  example 
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6.  VARIABLE  GLOSSARY 


A glossary  of  DYNASPHERE  variables  and  output  headings 
in  this  section.  Definitions  of  input  quantities  are  given 
description  section  of  this  report,  so  many  of  those  quanti 
found  here.  All  output  variables  and  headings  essential  to 
the  code  for  "production"  SGEMP  calculations  are  defined  in 


is  contained 
in  the  input 
ties  are  not 
the  usage  of 
the  glossary 


CLOSSAwr  OF  OVNISPMEHE  V*»UBLtS  AND  OUTPUT  HEADINGS 


ASPHER 

CURRENTS 

E 


EMISSION  AREA  (Mi) 

IN  2.D  PRINT,  actuallt  Should  be  current  oensitt 
ELECTRIC  FIELD  (VOLT/m) 


EMQ 

total  CHARGE  TO  BE  EMITTED  FROM  EACH  E“IS3I0N  ZONE  OVER 

THE  entire  emission  CURRENT  PULSE,  The  SUm  oe  These 

VALUES  Should  he  roughly  equal  to  the  integral  of  the 

EMISSION  CURRENT  density  OVER  time  4Nd  SPACE  FOR  The  ENTIRE 
PULSE,  AT  EACH  EMISSION  Time  step,  a particle  is  emitted 
FROM  Each  zone  with  an  AMOUNT  OF  CHARGE  EQUAL  TO 

TE  * £MQ  ID,  "HERE  TE  IS  DEFINED  BELON, 

energy 

EMISSION  electron  initial  energies  (REV), 

ER 

Radial  electric  field  (vplt/mi,  plotted  at  end  of  run 

FOR  VALUE  ON  inner  sphere  aT  TheTamO  VS,  Time. 

ER(0»R0oEG) 

ER  on  inner  sphere  RT  TheTaiO  and  RO  degree  positions 
(OR  closest  to  Them),  (vOlT/m) 

FGEnE 

qgene/oInjEk,  Fractional  amount  of  Emission  charge 
returning  TOhARO  sphere  at  each  radial  zone  (DIMENSIONLESS) 

H 

magnetic  field  (A“P/M) 

ITIME 

particle  time  step  number 

jnet 

approeimate  expression  for  The  net  Average  current  density 
leaving  the  inner  sphere,  equals  net  current  onto  inner 

SPHERE  DIVIDED  0Y  EMISSION  AREA,  AVERAGED  JNET  DESIGNATES 
jnet  also,  equals  JN0«  in  nun.SCL  cases,  (AmP/CMJ) 

JNOIk 

EMISSION  CURRENT  DENSITY  aT  THIS  TIME  AVERAGED  OVER 

EMISSION  AREA,  (AMP/CMJ), 

ma»,  no,  of 
particles  up  to 

MAXIMUM  VALUE  OF  NPART  UP  TO  TmIS  TIME, 

this  time 

NA 

number  of  angular  zones 

npacn 

CUMULATIVE  Number  of  particles  returning  to  inner  sphere 

UP  TO  this  time 
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NBtM> 

NLE*VE 

np*ht 

NR 

NRtT 

PHT 

POTf  N7UL 

0l 

Oi 

O3 

OBtC* 

8c.tNfc 

UHfRt 

aiN 

aiNji" 

QIN/OS 

UltiVE 

8NtT 

R 

RMiV 


BU4St»ST4TIC,  SINC.LE  SP-ERt  C 4LC11L4  T JCNS  CIM.*D  tOT4L  Nu“eiR 
CE  P4HIICLES  IN  S1UR46E  fl4N»  FOR  THOSE  P4RTICl.tS  OuTSlDE 
RESIGN  i«"E»E  GR£EN?S  FUNfTICN  FCR  » lELOS  IS  DEFINED. 

4 B4GI4L  FORCE  DUE  ’□  NET  CR4»&b  ON  3P"EHE  PULLS  ThesL 
P4RTItLES  BFCN  TO  Tm  REGION  R«E»E  ThE»  C4N  reenter  the 
SP4TI4L  “ESH  4N0  4G4IN  «4VE  bOT"  R4D14L  4ND  4»IiL  fields 
c‘lcul4ted  fbo“  The  greenis  function  4ct  on  the«*.  Those 
p»“TiClES  hITh  rinFTic  EnERc,»  gRFfteo  Ih4n  SPhERE  POTENTMl 
energy  tSC4PE  TO  infinity, 

CumI'LPTIvE  number  of  p4htici.es  rhICh  rave  cEft  system  up  to 
This  time  Hv  striking  inner  0«  cuter  80UN04RY 

Number  of  pahticles  being  follore.c  at  tris  tI"e.  includes 
Those  in  TrE  bank  in  OuaSI-STaTIC  calculations. 

Number  of  radial  cones 

8UASI.ST AT ic , single  SRrERE  CALCULATIONS  ONl,o  CUMULATIVE 
NUMBER  OF  particles  RETURNING  FROM  STORAGE  Bank  UP  TQ 
This  time  (SEE  NHANk) 

designates  INNER  SPHERE  POTENTIAL  RELATIVE  ’)  OUTER 
SPHERE  (VOLTS) 

electric  potential.  (YOlTS) 

generalized  COCIRPInate  CORRE  SRONCIng  to  radial  variable. 

generalized  coi'Rdinate  co““esronding  to  polar  angular 
variable  , 

generalized  coordinate  corresponding  T(i  4ZiruTh4u  angular 
variable. 

total  Particle  c-arge  ahich  hjs  siruc«  inner  sphere  ur 
to  This  time  (cnuL),  ALSU  used  TTi  designate  T-e  UGENE 
array,  oef ined  BELOr, 

jRRaY  nv  NR  values  giving  amount  of  charge  IN  EACH  Radial 
ZUNE  rhIch  is  returning  toraR'D  inner  sphere,  (COUL) 

charge  emitted  ERI'm  each  angular  Z NE  on  INNfR  sphere  minus 

charge  returned  to  each  anC.llAr 

ZONE  ON  Surface  of  inner  SRhERS,  return  C-argT  is  from 
RAhTICLES  I'Nly,  N!  ' -(E  T SiRFACF  Cl  rrENT  FlOr,  (touu) 

net  charge  out  In  A gIvEn  RADIUS  Including  t~e  net  ChAr&e  on 
the  inner  sphere  and  all  the  4>G,lAr  /ONES,  NOT  CALCULATED 
PROPERLY  AT  BRI5ENT,  (C'TIL), 

tctal  Charge  E“Itte,-  uR  tU  this  time  (coul) 

GIN/GNET,  GIVES  NET  Charge  out  TU  EACH  RADIAL  POSITION  AS 
fractuin  of  inner  sphere  net  Charge,  (DI“ensionless) 

NOT  calculated  properly  at  present, 

total  Charge  stri-inc.  spherical  surfaces  up  to  This  ti-e 
(CDUL),  includes  Charge  escaping  Past  Rmax  in  single 
sphere,  QUASIhSTaTIc  casts, 

N£ T Charge  on  IN'ER  sphere,  ESUALS  GI\JE«  • 8HAC«.  (COUL), 
ALSO  USED  AS  title  of  JheRE  ARRAY, 

Radial  zone  boundaries  OB  centers,  depending  on  KMICH  field, 
current,  or  charge  quantity  is  BEING  PRINTED  OUT, 
see  table  also  used  in  The  min, .PRINT  TO 

designate  Radial  r spions  of  the  'IRST  and  last  s particles 
presently  being  FOLlOhED,  also  uSf.o  in  emission  PaBTKLE 
characteristics  print  out  to  I'DICate  InPul  Radial 
pcspions  of  The  paryicles, 

cuter  RaOIaL  boundary  of  The  calculation,  can  be  a 
Radial  R iSPION  in  fmee  space  in  UuaSI-STatic  cases,  (m) 


ScUB 


SPEED 

8U»E*CE  C"*HGE 
OENsI T» 


SURF4CE  CD«BE^T 


SU«FCI'W(pS> 
90, US  DEg) 


T 


TE 


THEt» 

VB 

»The  T4 


SUBF4CE  CUBBtNTS  FLll»ISG  ON  INNjH  SPueBt  aT  f aCH  aNCUL** 

POSITION,  Ob’*1NED  fBO"  magnetic  field  JDST  outside 
SPhEBe  *hiCh  is  EV4LU4TEO  *T  EONE  CENTeBS  in  E*Cm 
COOHDIN4TE,  (AMPS) 

Emission  ElEctbqn  imtia^  speeds  (kev), 

SUBEacE  CmarGF  on  inner  and  UuTEb  spheres  1N  each  ZONE  ON 
THE  su^eacES,  inner  s‘’«f‘'E  ALftArs  pbINtEO,  OyfEB  sPnERE 
depends  on  I0PT(50),  The  production  print  flag,  obtained 
FROM  The  normal  ELECTRIC  EULD  NEAR  The  SURFACES,  ThE  SUR 
OF  The  values  on  The  inner  sphere  Should  agree  roughlt  hIth 
OnET,  the  net  Charge  on  the  inner  sphere  Obtained  from 
THE  particle  motion.  The  U1  VALUE  APPEARS  TO  ME  IN  ERROR 
FOR  This  particular  printout  OUANTITY,  (C'JUL), 

current  densities  flowing  on  inner  and  outer  spheres, 

INNER  SPHERE  ALMAYS  printed,  OUTER  SPHERE  DEPENDS  on 
icpt(3o),  The  production  print  flag.  These  current 
densities  are  obtained  from  The  magnetic  field  first 
ncnhZERU  value  near  The  spherical  surfaces.  The  value 
OF  Q1  appears  to  he  in  error  fob  This  particular  printout 
QUANTiTv,  (amp/m)  Plotted  nith  units  of  amps  at  end  of  run, 

surface  currents  on  inner  SP-ERE  closest  TP  the  u5,  90,  AND 
US  degree  positions,  RESPECTivEl*,  see  SCUB.  (Amps) 

The  f«AcT  angular  POSITIONS  CLOSES’  TO  The  “S,  90  AND  1 5S 

degree  locations  are  printed  out  in  The  Surface  current 
TI“F  history  Printout  near  the  last  page  UF  T"E  printout, 

TIME  (NSfC) 

fraction  of  THE  total  Emission  Current  pulse  area  uCcupied 
BY  THE  Current  Pulse  height  times  the  particle  time  s’er, 

(DIMENSIONLESS) 
polar  angle  (RADIANS) 
radial  velocity  (m/sEC) 
polar  velocity  (M/SEC) 
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